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ERRATA 


1*  PaSe  1?  Line  19  -  Read  "response"  instead  of  "respons" 

2.  Page  2,  Line  7  -  Delete  the  comma,  after  Z. 

3.  Page  4,  Line  lb  -  Nc;  c  is  lower  case. 

Pnge  5)  Line  12  -  Delete  comma  after  Cerenkov. 

5.  Page  6,  Line  8  -  Should  read  A  =r  99°N. 

6.  Page  8,  Line  22  -  Read  "three"  instead  of  "these". 

7.  Page  19}  Line  10  -  Insert  "is"  after  "this". 

8.  Page  19}  Line  19  -  Change  "or"  to  "of". 

9.  Page  23,  Line  13  -  Read  "energy"  instead  of  "enerby". 

10.  Page  29}  Line  2  -  Change  "alpha"  to  "alphas". 

11.  Page  30,  Line  7  -  Read  "monatomically"  instead  of 

'"monaomically".::-.  '  "  .  .  v- 

Page  3^}  Line  2b  -  Change  J(— E)  to  J(~E). 
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ABSTRACT 

The  flux  of  cosmic  ray  alpha  particles  has  been 
measured  at  ^  *  4l°N  and  55°N  and  the  energy  spectrum 
In  the  region  883-285  Mev/nucleon  has  been  determined  by 
the  use  of  a  combination  of  a  scintillation  counter  and 
Cerenkov  counter.  This  spectrum  Is  not  dependent  on 
geomagnetic  theory.  It  has  been  found  that  the  Integral 

i- 

alpha  spectrum  can  be  fitted  with  the  function 

J0  .  (E)  =  ^35 _  Particles  /M2-sec-ster 

*  (i+E)  ninrz 

and 

JQ  (E-i  )  =  88  iT  9  Particles  /M2-sec-ster  -  4l°N. 

J0^  (0.285  Bev)  =  292  ±  25  Particles  /M2-sec-ster  X=  55°N. 

J0^  (E)  Is  the  flux  of  ^  particles  at  the  top  of  the 
atmosphere  with  kinetic  energy  -  E  (in  Bev/Nucleon) . 

E^  is  the  cut  off  energy/nucleon  foroC  particles  at  4l°N. 

An  absorption  mean  free  path  =  45  £  7  g/cm2  has  been 
measured  for  alpha  particles.  The  properties  of  this 
new  detector  are  discussed  and  a  study  is  made  of  the 
energy  respons  of  the  Cerenkov  counter  to  fast 
multiply  charged  particles. 


* 
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Introduction 


Recently  it  has  been  demonstrated  1,2,3  that  it 


1.  J.  Llnsley ,  Phye.  Rev.  21>  1292  (1955). 

2.  N.  Horowitz,  Phys.  Rev.  £8,  165  (1955). 

3.  W.  R.  Webber  and  F.  B.  McDonald,  Phys,  Rev.  100 . 

1460  (1955). 


is  possible  to  obtain  consistent  and  accurate  flux 
values  of  the  low  Z,  components  of  the  primary  cosmic 
radiation  by  using  Cerenkov  counters  at  low  latitudes. 
These  experiments  have  provided  excellent  charge  reso¬ 
lution  for  Z  “  1,  2  and  adequate  resolution  in  the 
Z  “  3  -  8  range.  However,  it  appears  impossible  to  ex¬ 
tend  these  measurements  to  higher  latitudes  sihce,  with 
the  Cerenkov  counter,  slow,  high  Z  particles  give  pulses 
equivalent  to  those  of  fast  particles  of  lower  charge. 

If  such  an  extension  were  possible  it  would  be  of  great 
Interest  since  the  region  from  41°  to  55°  geomagnetic 
latitude  covers  the  velocity  sensitive  region  of  the 
average  Cerenkov  countermand  it  should  make  possible  the 
direct  determination  of  the  energy  spectrum  of  various 
components  without  having  the  uncertainties  introduced  by 
combining  the  geomagnetic  theroy  with  a  latitude  survey. 
In  addition,  an  electronic  detector  has  great  advantages 


Hill,  1949). 


It  will  be  noted  that  for  the  Cerenkov  detector  the 

( 

output  signal,  Nc ,  is  decreased  as  the  velocity  is  de¬ 
creased,  while  for  the  ionization  detector  the  output 
signal  is  increased  as  the  velocity  is  decreased.  In 
Figure  1,  a  plot  is  given  of  the  response  of  such  a  com¬ 
bination  of  these  two  types  of  detectors  to  protons 
and  alpha  particles  in  the  energy  range  above  320  Mev/ 
Nucleon.  The  ordinate  gives  the  most  probable  output 
signal  from  the  Cerenkov  counter  for  a  particle  of 
charge  Z  and  energy  E,  while  the  abscissa  gives  the 
corresponding  most  probable  energy  loss  <^>  in  the  ion¬ 
ization  device.  For  energies  less  than  1  Bev/Nucleon, 
the  curve  of  vs  Ne  can  be  represented  by:  k 

N0  -  Z2  Ki  -  K2  e 

where  Ki  and  K2  are  constants.  (3) 

With  the  individual  detectors,  resolution  of  the  cosmic 
ray  charge  spectrum  is  possible  only  for  relativistic 
particles.  However,  Fig.  1  demonstrates  that  the  com¬ 
bination  of  these  two  types  of  detectors  gives  charge 
resolution  which  is  independent  of  velocity.  This  re¬ 
solution  will  be  extended  below  the  Cerenkov  threshold 
due  to  material  in  the' telescope  beneath  the  ionization 
detector.  In  the  present  experiment  incident  particles 


must  have  an  energy  greater  than  120  Mev/Nucleon.  This 
prevents  slow  protons  from  extending  into  the  non- 
relatlvistic  alpha  region.  Since  the  response  of  both 
detectors  is  velocity  dependent,  it  follows  that  If 
effective  charge  resolution  can  be  achieved,  it  will  also 
be  possible  to  obtain  a  velocity  spectrum  of  the  various 
Z  components  of  the  cosmic  radiation  in  the  range  of  1  Bev 
to  100  Mev.  This  determination  of  the  energy  spectrum 
would  depend  only  on  the  properties  of  the  detector  and 
would  be  independent  of  geomagnetic  theory.  In  addition 
other  useful  Information  would  be  obtained.  Because  of 
the  directional  properties  of  the  Cerenkov^ counter  part¬ 
icles  lying  along  the  line  A-B  in  Fig.  1  would  represent 
the  upward  moving  albedo  of  the  “cosmic  radiation,  and 
it  would  be  possible  to  gain  some  idea  of  the  energy 
distribution  of  these  particles.  The  particles  l>ying  - 
along  R-B  where  R  is  at  the  vertical  geomagnetic  cut-off 
energy,  vould  represent  most  of  the  protons  in  the  re¬ 
turning  albedo  and  slow  secondaries  generated  In  the 
atmosphere.  r 

Detectors  using  this  principle  have  been  developed 
at  the  State  University  of  Iowa,  and  three  /lights  using 
Skyhook  balloons  have  been  completed. ■  Two  flights  were 
made  at  ^=4l°N  (San  Angelo,  Texas)  to  test  the  response 
to  fast  particles  and  to  note  whether  there  was  any 


background  in  the  multiply?  charged,  low  energy  particle 
region.  It  was  observed  that  while  the  charge  resolution 
was  comparable  to,  or  better  than,  that  obtained  by 
using  the  Cerenkov  counter  at  the  same  latitude,  it 
still  did  not  represent  an  optimum  response.  The  re¬ 
solution  was  further  Improved  by  better  selection  of 
photomultipliers  and  by  reducing  the'  amount  of  material 
in  the  telescope.  One  flight  was  made  at  >-=  i5°M 
(Minneapolis,  Minn.)  on  7  July  1955*  These  Initial 
flights  were  made  to  measure  the  proton  and  alpha  particle 
flux.  Further  flights  are  planned  at  58°,  55°  and  52° 
to  study  Z  =  1  -  6  components.  It  is  the  results  of 
the  flight  atX~  55°N//dealing  with  alpha  particles^ in 
addition  to  related  data  from  the  X  = flights^  that 
will  be  discussed  in  this  paper.  The  data  on  primary 
protons  and  cosmic  ray  albedo  will  be  discussed  in  a 
future  paper. 


Apparatus 


A  schematic  drawing  of  the  telescope  Is  given  in 
Figure  2.  A  scintillation  counter  using  a  Na  -  I 
(Th.  act.)  crystal  was  used  for  the  ionization  detector. 
Lucite  was  selected  for  the  Cerenkov  radiation  because 
its  index  of  refraction  is  convenient  in  the  energy 
range  of  greater  than  320  Mev/nucleon.  The  bottom 
telescopic  element  is  a  15.24  cm.  x  14.7  cm.  array 
of  Geiger  counters.  The  counters  have  an  effective 
length  of  14.7  i:  °*3  cm..  The  wall  thickness  is  0.088 
mm.  copper.  The  telescope  geometric  factor  is  4.02-^; 
0.15  cm .2;- g/fjer  for  Flight  III  and  7'25±Q3  cm2-ster 

for  Flight  I. 

Fig.  3  gives  a  schematic  drawing  of  the  electronic 
circuits.  The  arrangement  is  a  modification  of  the 
design  used  by  Ney  and  Thon^  on  the  first  Skyhook 


6.  E.  P.  Ney  and  D.  M.  Thon,  Phys.  Rev".  _81,  1068  (1951)- 

scintillation  counter  flight.  The  signals  from  the  two 
detectors  are  amplified  by  linear  amplifiers  and  dis¬ 
played  on  the  vertical  plates  of  two  cathode  ray  tubes. 

A  coincidence  is  formed  when  there  is  a  simultaneous 
pulse  from  the  Geiger  counter  tray  and  a  pulse  greater 
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than  .5  x  the  most  probable  energy  loss  for  a  proton 
at  minimum  ionization  in  the  crystal  channel.  The 
coincidence  pulse  operates  an  intensifier ' circuit  and 
the  pulses  from  the  two  scopes  are  recorded  by  an 
open  shutter  camera  with  continuously  moving  film. 

It  will  be  noted  that  this  is  Just  the  familiar  tech¬ 
nique  of  using  the  cathode  ray  tubes  as  master  co¬ 
incidence  elements.  The  pulses  from  the  lower  geiger 
counter  tray  also  actuate  a  shower  circuit  c  ■>  that 
there  is  an  indication  on  the  horizontal  CRT  plates 
when  an  event  occurs  and  more  than  one  counter  in  the 
lower  tray  is  discharged.  A  ring  of  gtfard  counters 
is  placed  above  the  lower  tray  to  assist  in  the  de¬ 
tection  of  stars  and  side  showers  and.  to  provide  better 
definition  of  the  telescope  geometric  factor.  These 
counters  also  give  deflections  on  the  horizontal  plates 
when  any  one  of  the  four  is  discharged.  The  action 
of  the  multiple  particle  detectors  is  to  give  the 
pulses  an  inverted  "J"  appearance  and  these  pulses 
will  be  referred  to  as  J  pulses  in.  contradistinction 
to  the  normal  pulses  arising  from  the  passage  of  a 
single  particle  through  all  these  elements  of  the  tele¬ 
scope.  For  flights  the  complete  assembly  with  the 
necessary  battery  packs  is  placed  in  a  thin  aluminum 
pressurized  container.  Pressure-altitude  data  is  ob¬ 
tained  by  photographing  a  model  FA160  Wallace  and  Tiernan 
gauge  in  an  Independent  unit. 
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III 

Detector  Resolution 

The  sea  level  cosmic  ray  mu  mesons  are  convenient 
particles  for  checking  the  resolution  of  the  telescope. 

In  actual  practice,  particles  of  a  given  charge  will 
not  lie  along  one  of  the  lines  shown  in  Figure  1,  but 
will  be  distributed,  in  a  region  about  this  line  due 
to  various  types  of  fluctuations.  In  the  scintillation 
counter  the  principal  fluctuation  is  the  Williams- 
Landau  effect.  This  fluctuation  arises  when  fast  charged 
particles  traverse  a  thin  absorber  since  large  energy 
transfers  to  single  electrons  can  occasionally  occur. 

The  most  probable  energy  loss^  £.p^is  given  by  Formula  1 
and  there  is  a  straggling  of  the  individual  losses 
around  this  value.  The  fluctuations  in  energy  loss 
were  first  calculated  by  E.  J.  Williams7  and  later 
more  accurately  in  closed  form  by  Landau^  and  by  Symon^. 
The  curves  presented 


7.  E.  J.  Williams,  Proc.  Ruy,  Soc.  A12 5,  420  (1929). 

8.  L.  Landau,  Journal  of  Physics  USSR  8,  201,  (1944). 
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in  Symon’s  thesis  are  used  in  the  calculation  on  fluctu¬ 
ations  in  energy  loss.  Fig.  4  gives  an  observed  pulse 
height  distribution  for  sea  level  mu  mesons  and  the 
predicted  curve  based  on  the  Williams-Landau  effect. 

It  is  seen  that  excellent  agreement  with  the  predicted 
curve  is  obtained  on  the  high  energy  side'  of  the  curve. 
The  effect  of  electrons  escaping  from  the  crystal  is 
partially  offset  by  the  presence  of  1.5  g/cm  of  glass 
and  A1  above  the  crystal.  In  the  calculated  curve 
a  mean  meson  energy  of  1  Bev  is  assumed.  The  actual 
mu  meson  energy  spectrum  introduces  a  negligible  error 
in  the  calculation.  The  additional  broadening  of  the 
distribution  is  due  primarily  to  non-unif ormity  of 
light  collection  with  a  much  smaller  contribution  from 
photon  statistics.  In  Figure  5  is  shown  the  calculated 
Williams-Landau  distribution  for  o  particles  with 
energies  of  4,  1.5?  .8,  .4,  .3  and  .2  Bev/nucleon.  The 
distributions  become  more  gaussian  in  shape  and  become 
progressively  sharper  as  the  kinetic  energy  of  the 
particles  is  decreased.  A  fortunate  circumstance  oc¬ 
curs  in  that  the  fast  a  particles  (E  >  1.5  Bev.  nucleon) 
and  sea  level  mu  mesons  will  both  have  approximately 
similarly  shaped  distributions  in  the  ionization  de¬ 
tector.  This  arises  because  the  two  major  sources  of 
fluctuation  -  the  Williams-Landau  effect  and  the  non¬ 
uniformity  of  light  collection  -  are  the  same  for  both. 


The  fractional  spread  of  the  Wllliams-Landau  distri¬ 
bution  is  a  very  slowly  varying  function  of  £ for 
fast  particles  when  E/zm/’’  >  2.  and  E  >>  £  , 

p 

A  comparison  of  the  sea  level  mu  mesons  distribution, 
shifted  to  the  -region,  and  primary  cosmic  ray  primary 
alpha  particle  distribution  obtained  on  Flight  #1, 

January  1955  at  X  =  4l°  is  given  in  Figure  6.  As  will 
be  apparent  later,  the  major  correction  for  the  raw 
data  involves  the  high  energy  tail  on  the  fast  particle 
distributions.  The  close  agreement  between  the  predicted 
distribution  and  the  experimental  distribution  as  shown 
in  Fig.  4  and  6  give* increased  confidence  that  this 
correction  can  be  applied  with  a  high  degree  of  ac¬ 
curacy. 

In  the  Cerenkov  counter  the  major  source  of  fluc¬ 
tuations  is  due  to  photon  statistics.  While  the  number 
of  useable  photons  produced  by  the  passage  of  a  fast 
singly  charged  particle  through  the  Incite  is  of  the 
order  of  120Q,  the  photoelectric  efficiency  of  the 
photocathode  surface  is  low  and  the  actual  number  of 
photoelectrons  will  be  of  the  order  of  40  resulting  in 
an  approximate  Poisson  distribution.  Fig.  7  gives  the 
experimental  Cerenkov  distribution  for  sea  level  mu 
mesons.  The  dotted  line  represents  the  Poisson  dis¬ 
tribution  for  n  =  34.  Since  part  of  the  broadening 


of  the  distribution  is  due  to  cascading  stages,  each 

o 

having  a  poisson  distribution,  the  actual  number  of 
photoelectrons  -will  be  greater  than  H  =  3^  and  will 
depend  on  the  secondary  emission  ration  of  the  various 
stages  and  focusing  losses  between  dyno.de e.  It  is 
seen  in  Figure  6  that  the  observed  distribution  is 
■■asymmetric  .  This  represents  an  effect  in  the 
Cerenkov  counter  similar  to  the  Willlams-Landau  effect 
in  an  ionization  detector. 
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IV 

FLIGHT  DATA 

The  time-altitude  data  for  Flight  I  at  X  =  4-1°  (San 
Angelo,  Texas)  on  17  Jan.  1955  and  Flight  III  at  X  =  55° 
(Minneapolis,  Minn.)  on  7  July  1955  are  given  in  Figure  8. 
The  time  in  each  case  is  central  standard  time.  The 
accuracy  of  pressure  determination  is  better  than  . 5mb. 

The  agreement  with  the  telemetered  data  furnished  by  the 
balloon  contractor,  Winzen  Research,  Inc.  of  Minneapolis, 
Minnesota  was  generally  within  1  mb. 

In  both  flights  the  equipment  was  recovered  in 
operating  condition  and  immediate  tests  were  carried 
out  to  insure  that  the  apparatus  was  functioning  properly 
and  in  particular  that  the  geiger  counter  plateaus  had 
not  shifted.  In  each  case  there  was  no  discernible 
shift  in  any  component. 

Reeuction  of  the  data  was  carried  out  by  displaying 
the  film  on  a  Model  C  Recordax  microfilm  reader  and, 
for  each  event,  measuring  the  pulse  height  in  the  Cerenkov 
channel  and  in  the  ionization  channel.  It  was  also 
noted  whether  one  of  the  multiple  event  detectors  had 
been  actuated  for  a  given  event.  Appropriate  time 
intervals  were  selected,  with  a  maximum  interval  of 
20  minutes  being  used.  The  data  were  tabulated  on  a 
4-0  x  40  rectangular  array.  On  the  film,  the  average 
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spacing  between  counts  is  not  large  and  care  is  required 
to  insure  that  the  proper  pair  of  pulses  is  selected. 

In  less  than  1%  of  the  counts,  it  was  not  possible  to 
make  an  effective  separation  between  a  pair  of  events. 

After  the  film  data  had  been  transferred  to  the 
40  x  40  array,  it  became  obvious  that  due  to  the  widths 
of  the  distributions,  it  would  not  be  possible  to  use 
equation  3  in  obtaining  the  flux  values.  It  was  then 
decided  to  divide  the  data  into  various  ionization 
intervals  and  study  the  Cerenkov  distributions  in  these 
intervals.  This  is  illustrated  in  the  corrected  data 
of  the  particle  region  given  in  Fig.  9.  These 
give  the  distribution  of  the  regular  Cervenkov  counts 
in  the  selected  ionization  intervals.  The  same  data, 
along  with  the  distribution  of  the  "J"  counts  are  given 
in  Fig.  10.  No  corrections  have  been  applied.  The 
ionization  intervals  chosen  correspond  to  kinetic 
energies  per  nucleon  at  the  top  of  the  atmosphere  of 
>1.5  Bev.,>.883  Bev.,  -.570-. 883  Bev. ,  .436-  .630  Bev., 
.320  -  .436  and  .286  -  .361  Bev..  The  intervals  were 
selected  so  that  there  were,  to  a  crude  approximation, 
equal  numbers  of  alpha  particles  in  each  of  the  inter¬ 
vals  which  were  in  the  energy  sensitive  region  (i.e., 

E  900  Mev/nuc).  Certain  features  of  the  curves 
are  immediately  apparent,.  One  notes  in  the  fast 
particle  region  the  excellent  peak  to  valley  ratio, 
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and  the  manner  in  which  the  curves  go  to  zero  on  the 
high  side  of  the  distribution.  It  is  also  seen  in 
Fig.  10  that  the  J  distribution  is  slightly  peaked  under 
the  alpha  region  and  that  most  of  these  counts  can  be 
accounted  for  by  6  rays  and  interactions  in  material 
in  the  telescope  below  the  Cerenkov  detector.  It  is 
felt  that  the  good  peak  to  valley  ratio,  the  sharp 
fall-off  on  the  high  side  of  the  distribution  and  the 
MJH  pulse  distribution  indicate  a  minimum  of  background 
in  the  alpha  region.  This  further  strengthened  by  the 
Observed  altitude  dependence  given  in  Fig.  11  which 
will  be  discussed  later. 
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CORRECTIONS  IN  ALPHA  REGION 

There  are  four  corrections  to  be  applied  to  the 
observed  alpha  data  in  order  to  deduce  an  absolute  flux 
at  the  top  of  the  atmosphere.  There  are l  a  general 
background  correction,  a  correction  for  the  absorption 
of  alpha  particles  in  the  telescope  and  in  the  atmos¬ 
phere  above  the  telescope  due  to  nuclear  interactions, 
a  £  ray  correction,  and  a  correction  for  the  production 
of  alphas  in  the  atmosphere  by  the  fragmentation  of 
primary  heavy  nuclei. 

The  shape  of  the  background  correction  is  taken 
to  be  the  interval  at  3  x£p  for  protons  at  minimum 
ionization.  In  Fig.  1  this  corresponds  to  a  line 
drawn  parallel  to  the  ordinate  and  Intersecting  the 
abscissa  at  3*6  on  the  arbitrary  scale.  The  ionization 
distribution  rises  rapidly  so  it  is  possible  to  select 
a  region  Just  behind  a  group  of  particles  with  the  de¬ 
sired  charge  and  this  region  will  be  almost  completely 
free  of  particles  with  this  charge.  The  flight  at 
X=4l°N  demonstrates  that  the  shape  of  the  background 
curve  did  not  change  across  the  complete  alpha  interval. 
The  background  curve  obtained  at'X=  55°N  was  normalized 
on  the  low  energy  side  of  the  Cerenkov  distributions. 

The  corrections  to  the  various  intervals  by  this  method 
are  given  in  Fig.  10. 
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In  order  to  correct  for  the  attenuation  of  the 
alpha  particles  by  nuclear  interactions  it  is  nec¬ 
essary  to  know  the  interaction  mean  free  path,  . 

There  are  several  experimental  determinations  of  . 

Waddington10  has  used  emulsion  at  ^  =  55° »  Webber 

10.  C.  J.  Waddington,  Phil  Mag.  4£,  1312  (195*0. 


and  McDonald3  at  ^  =  41°  observed  the  atmospheric 
attenuation  as  detected  by  a  Cerenkov  counter  and 
Quarentl  and  Zornll  used  90  Mev/  nucleon  alpha 


11.  G.  T.  Quarentl  and  G.  T.  Zorn,  Nuovo  Cimento  1, 
1282  (1955). 


particles  from  the  Berkley  cyclotron.  All  values 
obtained  are  in  good  agreement  with  the  semi-empirical 
relation  proposed  by  Bradt  and  Peters-*^: 


$ 


12. 


Bradt  and  Peters  Phys.  Rev.  77,  5*^  (1950). 

<T  -  7T(r,  +  &  -  £*r)e 


£  =  /.45  x/0-,3A^  cm 
a r  -  ,8S  * /0~/s  cm. 


A]_  =  mass  number  of  Incident  nuclei 
A£.  -  mass  number  of  target  nuclei 
-  interaction  cross  section 


This  gives  a  value  of  ^^(air)  =  44  g/om2.  The  observed 
atmospheric  absorption  of  Of  's  in  this  experiment  is 
given  in  Fig.  11.  This  gives  a  value  of 
^  =45±7g/cm2.  The  work  of  Qwarenti  and  Zorn  at  low 
energies  gives  added  confidence  that  the  collision  cross 
section  is  not  strongly  energy  dependent.  A  compilation 
of  the  observed  values  of  Aq<  for  the  various  experiments 
is  given  in  Table  I.  The  nuclear  interaction  correction 
factor  for  Flights  I  and  III  is  summed  up  in  Table  II. 

The  (fray  correction  is  necessary  as  it  is  possible 
for  particles  traversing  the  top  two  elements  of  the 
telescope  to  produce  knock-on  electrons  which  will 
trigger  one  of  the  gelger  counters.  For  particles  which 
would  normally  pass  through  the  telescope,  this  results 
in  their  being  classified  as  J  counts.  Some  particles 
which  miss  the  bottom  element  will  appear  as  regular 
counts  and  there  will  be  a  spurious  increase  of  the 
geometric  factor.  Because  of  the  straggling  of  the 
electrons  in  the  block  and  in  the  counter  walls,  it  is 
difficult  to  compute  this  correction  with  confidence. 
However  an  upper  limit  to  the  correction  can  be  obtained 
and  this  upper  limit  probably  represents  a  good  ap¬ 
proximation  to  the  true  correction.  If  one  neglects 
straggling  then  the  number  of  knock-one,  per 

particle  which  trigger  a  counter  in  the  lower  tray 


Table  I 


Observed  M.F.P.  of  Helium  Nuclei  in  Emulsion  and  Air 


Investigator 

Lat.° 

Detector 

(4?) 

g/ cm^- 

from  eq 
g/cim 

Waddington10 

550 

Emulsion 

6o±l7 

72  . 

Quarenti 

Zorn11 

90  Mev/Nuc 
©(from  Berkeley 
cyclotron 

Emulsion 

80  ±6 

72 

Webber 

McDonald3 

410 

Cerenkov 

counter 

43  X  8 

44  • 

McDonald 
(Present 
Experiment ) 

55° 

Cerenkov 
and  sclnt. 
Counter 

45-7 

44 

rv>* 


Table  II 


Flight 

I 

III 


Correction  for  Attenuation  ofoc's 
by  Nuclear  Interactions 


Lat. 
4 1°N 
55°N 


*1 

Atmos. 

Depth 

of 

Detector 


*2 


Amount  of 
Material  in 
Te le scope 

2r 


e-teL  +  \ 


15.1  g/ cm2  9.08  g/cm2  .580 


12.4  9.2  -813 


can  be  calculated.  This  quantity  can  also  be  deter¬ 
mined  experimentally  for  sea  level  mu-mesons,  and 
this  is  designated  as  (N^)  exp.  The  ratio  of 

exp/  (N<}x/)c  =  Y  is  taken  as  a  correction  factor 
to  be  applied  to  (%t><)c.  If  one  neglects  straggling 
then  the  number  of  knock-ons  which  trigger  the  bottom 
tray  is  given  by 


)#<* 


dds 


JK,  SsvtW  ls  (5) 

The  range-energy  relation  for  electrons  is  given  by^-3 
E  =  ax  +  c  =  1.83x  +  .133  Mev. 

E-^  =  ax^  +  c  =  1.25  Mev. 


and 


13.  Glendenin,  L. ,  Nucleonics  2,  No.  1,  12  (194-8). 

In  these  formulas  p  and  Z  are  the  velocity  and  charge 
of  the  incident  particle,  x-^  is  the  thickness  or  the  geiger 
counter  wall,  X2  is  the  total  amount  of  material  in  the 
telescope  and  E2  is  equal  to  15  Mev.  E g  max  is  the 
maximum  energy  which  can  be  transferred  to  an  electron 

by  a  heavy  particle  of  velocity  p. 

/  c 

If  F^  max  =  x2  +  C  =  E2  (this  corresponds  to 

3  Bev/nuc  for  alpha  particles)  then  (N.  )  reduces 

*  z  c 

to 


€Zz 


v  csmdK/  <-/  C6/7M  _ 
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r*  ^ 

The  experimentally  determined  ratio\  Y*  was  used,  to 
correct  (  ^  ^  )  at  various  energies.  These  cor¬ 
rections  a.re  given  in  Table  III.  It  is  seen  that  in 
the  energy  sensitive  region  of  the  detector  (l.e.,  H" 
1  Bev),  this  correction  is  negligible.  For  the 
Flights  at  41°  the  bottom  ring  of  guard  counters 
was  omitted  and  the^-ray  correction  was  -2o/o. 

In  this  correction  the  $  -ray  angular  distribution 
of  cos^ ^©measured  by  Brown,  McKay  and  Palmatier^ 


14.  W.  W.  Brown,  A.  S.  McKay  and  E.  D.  Palmatier, 
Phys.  Rev.  £6,  50 6  (1949). 


on  sea  level  mesons  was  used,  where  B  is  the  angle 
between  the & -ray  emerging  from  the  block  and  the 
track  of  the  particle  which  produced  the  -ray. 

To  extrapolate  the  observed  alpha  flux  to  the 
top  of  the  atmosphere  it  is  also  necessary  to  correct 
for  the  production  of  energetic  e<  particles  by  the 
fragmentation  of  heavy  nuclei  above  the  telescope. 


Table  III 


E_ 

Kinetic  Energy/ 
Nucleon  of 
Incident  o( 

/ 

Er 

o  max 

Maximum  Energy 
$  Ray  which  can 
be  produced  by 
alpha  with  energy 

Number  of 
^Rays/oC 
which  trigger 
g-elger  counter 
array 

(Straggling 

neglected) 

=E 

Y 

number  of 
(fraye/of  which 
trigger  geiger 
counter  array 
with  empirical 
straggling 
correcti on 

4  Bev 

28  Mev 

.18 

.045 

2 

8.75 

.10 

.025 

1 

3.29 

.037 

.010 

.8 

2.4  2 

.017 

.043 

.  6 

1.7 

.005 

.001 
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This  Is  similar  to  the  correction  necessary  in  the  Li, 

Be  and  B  region.  Noon  and  Kaplonl5  have  derived  the 

15.  J.  H.  Noon  and  M.  F.  Kaplon,  Phys.  Rev.  97,  lAl  (1955). 


fragmentation  probabilities  for  incident  medium  and 
heavy  nuclei  in  air  based  on  the  observed  Interactions 
in  emulsion.  The  value  for  the  fragmentation  prob¬ 
ability,  P.,  ^  which  is  the  average  number  of  alphas 
produced  in  a  nuclear  interaction  initiated  by  the  xth 
component  and  appropriate  heavy  flux  values  are  summed 
up  in  Table  IV.  A  standard  diffusion  equation  Is  used 
to  obtain  the  actual  number  of  secondary  alphas  Incident 
on  the  telescope.  t 
* 

At  A  *  55°  number  of  secondary  alphas  at  12  g/cm2  = 

11  %  4  particles/m^-sec-ster . 

At  X=  4l°  number  of  secondary  alphas  at  15  g/cm^  = 

6.5  i  3  particle s/m-2-sec-eter. 

This  correction  has  not  been  applied  to  the  differential 
and  integral  spectrum  in  Fig.  14  and  15  and  in  Table  4. 


Plux  Values  and  Fragmentation  Probabilities  Used  in  Fragmentation  Correction 
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V 

ENERGY  CALIBRATION  ■ 

The  selected  intervale  of  ionization  were  cali¬ 
brated  in  terms  of  the  most  probable  energy  loss  for 
particles  at  minimum  ionization,  (6^)^  min' 

From  the  data  at  41°  it  is  known  that  some  25% 
of  the  alphas  at  55°  will  have  energies  greater  than 
2  Bev/Nuc  and  will  be  at  minimum  ionization.  The 
ionization  distribution  of  these  particles  can  be 
studied  by  selecting  an  Interval  along  the  high  side 
of  the  Cerenkov  distribution.  In  Fig.l  this  corre¬ 
sponds  to  taking  an  interval  parallel  to  the  £p  axis 
and  intersecting  the  Cerenkov  axis  at  21.  This  gives 
a  pulse  height  distribution  similar  to  that  shown  in 
Fig.  6.  The  peak  in  the  distribution  thus  obtained  is 
(Ep)g<  min  and  the  values  of  (Ep)c<  of  the  other 

intervals  are  known  in  terms  of  (Ep)^mj_n  since  the 
system  is  a  linear  one.  The  peak  in  the  ionization 
distribution  is  quite  accurately  defined  and  the 
position  of  (£p^min  the  40  x  40  array  is  known 

to  3%.  A  curve  of  (£p)o<  vs.  kinetic  energy  as  cal¬ 
culated  by  Eq.  1  is  given  in  Fig.  12.  The  abscissa 
could  also  have  been  plotted  in  terms  of  (  £p)<A  / 

(E.p^  min*  si-nce  this  was  the  quantity  actually  used 
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VI 

ALPHA  ENERGY  SPECTRUM 

To  derive  an  energy  spectrum  from  the  distributions 
given  in  Fig.  10  it  is  necessary  to  relate  the  observed 
number  of  particles  within  an  Interval  of  energy  loss 
to  the  actual  number  in  the  corresponding  energy  inter¬ 
val.  The  calibration  of  the  energy  intervals  in  terms 
of  ionization  loss  was  discussed  in  the  previous  section. 
However,  due  to  fluctuations  introduced  by  detectors, 
the  number  of  particles  within  an  interval  of  energy 
loss  is  not  necessarily  the  number  in  the  related  en¬ 
ergy  Interval.  The  largest  correction  will  be  the 
Wllliams-Landau  correction  for  high  energy  alphas. 

It  is  a  major  correction  only  in  the  energy  region 
1=  1  Bev  and  becomes  completely  negligible  at  energies 
less  than  .6  Bev./nuc.  This  is  illustrated  in  Fig.  5. 

As  has  been  discussed  in  Section  3  this  correction 
can  be  computed  quite  accurately.  This  correction  has 
been  calculated  for  intervals  I  and  II  (  E  >  883  Mev/ 
nuc)  and  III  (570-883  Mev/nuc.).  The  effect  of  ap¬ 
plying  this  correction  is  to  remove  the  asymmetry  in 
the  ionization  distribution.  Thus  the  pulse  height 
distributions  from  both  detectors  should  be  approxi¬ 
mately  gaussian  in  the  energy  sensitive  region.  If 
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now  the  Intervals  have  been  selected  so  that  there  are 
equal  numbers  of  alpha  per  interval  after  the  Williams- 
Landau  asymmetry  has  been  removed,  then  this  number 
should  represent  a  good  approximation  to  the  actual 
number  of  particles  -in  the  energy  interval  since  the 
loss  to  neighboring  intervals  will  tend  to  be  balanced 
by  the  gain  from  these  intervals.  The  energy  spectrum 
data  are  summarized  in  Table  V.  It  is  seen  that  this 
condition  has  been  well  satisfied  in  intervals  3  and  k 
but  the  value  of  falls  off  in  intervals  5  and  6. 

A  correction  factor  was  computed  by  the  synthesis  of 
the  pulse  height  distributions.  This  correction  factor 
was  never  greater  than  10# .  This  gives  added  confidence 
that  this  method  of  analysis  has  not  Introduced  major 
errors  in  the  derivation  of  the  energy  spectrum.  In 
a  similar  manner  the  contributions  of  intervals  I  and 
II  to  III  were  computed.  This  is  necessary  because  of 
the  great  difference  in  in  the  3  intervals. 

These  corrections  have  all  been  applied  with  the 
Williams -Landau  correction  in  column^  of  Table  V. 

In  Column  7,  No<  is  converted  to  a  vertical  flux  value 
and  in  Column  8  the  <$"  ray  correction  and  the  correction 
for  nuclear  absorption  in  the  telescope  and  in  the  at¬ 
mosphere  above  the  telescope  have  been  applied.  In 
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Column  9  the  values  for  the  differential  energy  spec¬ 
trum  are  given.  The  flux  data  for  Flight  I  ,are  also 
included  in  Table  V.  Fig.  13  and  Fig.  14  give  the 
points  obtained  for  the  differential  and  integral 
energy  spectrum.  In  each  graph  the  high  energy  point 
makes  use  of  at  4l°N.  A  geomagnetic  cut-off  of 

1.65  Bev/Nuc  has  been  calculated  for  this  point.  This 
is  the  only  point  dependent  on  geomagnetic  theory. 

Also  in  Fig.  13  are  the  experimental  points  measured  by 
Waddington10  at  A  =  55°  using  emulsion.  The 
general  agreement  is  good.  The  data  in  Fig  13  and  14 
have  been  fitted  with  a  curve  of  the  form 

J0fl(  (£e)  =  (1+E7Y 

E  =  kinetic  energy  in  Bev. 

It  is  found  that  a  best  fit  is  found  for: 

y  =  1.4 +.2  and  c  =  431. 

The  vertical^  flux  values -at  A  =  are 

J0  (E>285  Mev/Nuc )  =  3°6~24  particular/M2-sec~ster  (without 
°<  fragmentatlor 

ft  ft  //  correction) 

Jo*  (El)  =  96±  9 

Ex  =  vertical  geomagnetic  cut  off  in  Bev/Nuc  at  41°. 

When  the  fragmentation  correction  is  applied 

JQ_  (E  285  Mev/Nuc)  =  292+  25  Particles/M2-sec-ster 

l )  /)  n  1 1  H 

Jo*  (El)  =88+9 
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Because  of  a  gain  shift  caused  by  repairs  Just 
prior  to  launching,  it  was  not  possible  to  say  anything 
concerning  the  geomagnetic  cut-off  at  ^  =  55°  except 
that  it  is  less  than  285  Mev/Nuc. 

Table  VI  provides  a  summary  of  o<  flux  measurements 
at  ^  =  41°  and  55°»  The  over-all  agreement  is  quite 
good. 


16.  Perlow,  Davis,  Kissinger  and  Shipman 
Phys.  Rev.  88,  321  (1952). 

17.  J.  Linsley,  Phys.  Rev.  101,  826,  (1956) 

18.  L.  Bohl,  University  of  Minnesota,  Thesis  (195^)* 

19.  L.  R.  Davis,  H.  M.  Caulk  and  C.  Y.  Johnson 
Phys.  Rev.  101,  800  (1956). 

20.  F.  B.  McDonald,  University  of  Minnesota  Thesis 

(1955). 
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Climax,  Colorado.  There  are  no  observed  fluctuations 
during  the  period  of  the  flight  and  the  neutron  counting 
rate  is  in  close  agreement  with  the  yearly  mean  -  illus¬ 
trating  that  this  was  a  “normal  cosmic  ray  day". 
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VIII 

ENERGY  RESPONSE  OF  CERENKOV  COUNTER  TO  FAST 
MULTIPLY  CHARGED  PARTICLES 

The  theoretical  Z  and  J3  dependence  of  the  Cerenkov 
radiation  is  given  hy  equation  2.  However  this  equation 
was  not  used  explicitly  in  deducing  the  energy  spectrum. 
It  was  assumed  merely  that  the  Nc  was  Z  dependent  and 
decreased  monaomically  as  the  velocity  decreased.  It 
is  of  interest  to  measure  the  Z  and  f}  response  of  the 
Cerenkov  counter  to  fast  multuply  charged  particles. 

By  comparing  the  relativistic  proton  and  alpha  dis¬ 
tribution  (E  2  Bev/Nuc)  it  was  established  that  for 
/ 3^1,  the  radiation  varied  as  Z2  within  the  7%  ex¬ 

perimental  error.  The  Q  dependence  of  the  Cerenkov 
radiation  can  be  studied  by  noting  the  position  of  the 
peak  in  the  Cerenkov  distribution  in  the  various  ioni¬ 
zation  intervals  given  in  Fig.  10.  These  experimental 
points  have  been  normalized  to  the  theoretical  curve 
at  interval  I  (E  =  3  Bev.).  Agreement  with  the 
theoretical  curve  is  excellent  for  the  higher  values 
of  E.  However  at  levels  near  the  Cerenkov  threshold 
there  appears  to  be  some  fluorescence.  This  is  not 
unexpected. 
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CONCLUSIONS 


It  is  believed  that  the  method  described  in  this 
paper  is  equally  applicable  to  the  proton,  light  and 
medium  cosmic  ray  components.  It  is  interesting  to 
note  that  at  300  Mev/Nuc,  the  alpha  spectrum  is  still 
rising  rapidly  toward  lower  energies.  It  is  planned 
to  extend  these  measurements  to  higher  latitudes'-- in 
the  near  future. 
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Fig. 


Fig.  2 

Fig.  3 

Fig. 


Fig.  5 


Fig.  6 


Fig. 


FIGURE  CAPTIONS 


.  Plot  of  Cerenkov  Pulse  Height  in  1"  Lucite 
radiator  vs  the  most  probable  energy  loss  for 
protons  and  alphas  in  a  1/8”  Na-I  crystal. 

Fast  albedo  particles  ■will  be  along  A-B. 

R-B  represnets  particles  below  geomagnetic  cut¬ 
off.  Triangles  represent  observed  response  of 
Cerenkov  detector  to  fast  multiply  charged 
particles  in  Flight  III. 

.  Schematic  Drawing  of  Telescope. 

.  Schematic  Diagram  of  electronic  circuits. 

.  -  0  -'Experimental  pulse  height  distribution  for 
sea  level  mu  mesons.  Dashed  line  is  calculated 
curve  for  Williams-Landau  distribution. 

.  Calculated  Williams-Landau  Distribution  for  alpha 
particles  of  1.5,  -8,  .4,  .3  and  .2  Bev  Nucleons. 
P  (E)  represents  probability  of  energy  loss  E  in 

crystal  by  a  particle  of  given  energy. 

/  , 

.  Comparison  of  alpha  pul/se  height  distribution 
(-  0  -)  in  ionization  detector  (Fit.  I,  \  -  41°) 
with  sea  level  mu  meson  distribution.  (Both 
curves  normalized  to  same  area). 

.  Experimental  sea  level  Mu  meson  pulse  height 
distribution  in  Cerenkov  detector.  Dashed  line 
represents  Poisson  distribution  with  n  =  34. 
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Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


Fig. 


8.  Time  altitude  curve  for  Fit.  I  (San  Angelo, 

Texas,  17  January  1955)  and  Fit.  Ill  (Minne¬ 
apolis,  Minn.,  7  July  1955)* 

9.  Raw  Differential  Pulse  Height  Data  in c<  region. 

Ionization  intervals  have  been  selected  so  that 
there  is,  to  a  crude  approximation,  equal 
number  of  o(  counts/interval  in  the  energy 
sensitive  region.  The  ordinate  scale  represents 
the  actual  numbers  of  counts  at  altitude. 

10.  Differential  pulse  height  data  of  Figure  9* 

Heavy  dashed  line  represented  background  counts 

in  interval.  - & - gives  distribution  of  J 

counts  in  interval. 

11.  Vertical  intensity  of  alpha  particles  as  a 

function  of  atmospheric  depth.  This  gives  an 

*» 

apparent  mean  free  path  of  45  ~  7  g®s  as  in¬ 
dicated  by. solid  line. 

12.  Plot  of  Kinetic  Energy/Nucleon  a.s  a  function 

of  the  most  probable  energy  loss  in  the  crystal 
for  of  particles . 

13.  Differential  Energy  spectrum  for  primary  alphas. 

The  emulsion  data  obtained  by  Waddington^  is 
shown  for  comparison.  The  data  has  been  fitted 
to  a  curve  of  the  form  J  (  =  E)  =  g)Y  • 

The  point  on  the  right  at  1.3  Bev  was  obtained 
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using  data  at  ^  -  41°.  This  is  only  point 
involving  geomagnetic  theory  and  was  disre¬ 
garded  in  fitting  the  data. 

Fig.  14.  Integral  energy  spectrum  of  primary  alpha  part¬ 
icles.  Point  at  1.7  Bev/  nucleon  based  on 
\  =  41°  data. 

Fig.  15-  Observed  hourly  rates  for  total  alpha,  fast 

alpha  (E  >  800  Mev/nucleon)  slow  alpha  (E<  800 
Mev/Nucleon)  and  total  telescope  rate.  Time 
is  Central  Standard  Time. 

Fig.  16.  Plot  of  Cerenkov  light  in  kinetic  energy/ 
nucleon  for  alpha  particles.  Experimental 
points  have  been  normalized  to  theoretical 
curve  at  highest  experimental  point. 
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